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Abstract 



We have studied the exciton states in two coupled conjugated polymer chains 
which are modeled individually by the Su-Schrieffer-Heeger Hamiltonian and 
coupled by an interchain electron-transfer term. Both the intra- and inter- 
chain long range Coulomb interactions are taken into account. The proper- 
ties of the lowest symmetric and antisymmetric exciton states are extensively 
discussed for both the parallel and anti-parallel ordering between these two 
chains. It is found that, for these two kinds of ordering, the features of excitons 
are quite different. Possible implications for the experiment of luminescent 
polymers are also addressed. 
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I. INTRODUCTION 



The significance of excitons in conjugated polymers lias been recognized for many years 
in a class known as polydiacetylene Recently, the interest in excitons has been height- 
ened by discovering that poly(]9-phenylene vinylene) (PPV) and its derivatives can be used 
as the active luminescent layer in electroluminescent light-emitting diode devices [0, since 
it is believed that radiative recombination of singlet excitons gives rise to luminescence. 
Although a general picture of understanding of photoinduced absorption (PA) and pho- 
toconductivity (PC) experiments of PPV and its derivatives remains a subject of intense 
debate [0HTT||, the dramatically different PA behavior of dilute solutions and thin films of 
poly[2-methoxy, 5-(2' ethyl-hexoxy)-l,4 phenylene vinylene] (MEH-PPV) [|^, together with 
the photoconductivity being observed to start at the absorption edge in PPV ||13|, clearly 
indicate the important role the interchain coupling plays. Actually, electron-diffraction ex- 
periments indicate that the interchain coherence length of PPV is about 60 A [Q. However, 
the properties of excitons in coupled chains, to the best of our knowledge, have never been 
discussed. 

There have been several works to develop the exciton theory in a single polymer chain 



||r5Hl8[ . Abe and co-workers introduced the standard exciton theory |T9| to one-dimensional 
polymers |T^. Within the single configuration interaction (SCI) approximation, the energy 
levels and wave functions of exciton states in a long chain can be obtained. In this paper, 
we will extend their approach and explore the features of excitons in two coupled polymer 
chains. Interestingly enough, we find that, for two kinds of ordering, i.e., the parallel- 
and anti-parallel- ordering between the two chains, the properties of excitons are quite 
different. Although this study does not result in quantitative explanations of the PA and 
PC experiments, it is helpful for understanding the photophysics better in luminescent 
polymers like PPV, and is also a foundation of further calculations of optical properties of 
coupled polymers. In Sec. II, the model is defined and the formulation is derived. In Sec. 
Ill, we present numerical results on the lowest symmetric and antisymmetric excitons and 
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discuss some implications for experiments. 



II. FORMALISM 



We start with the two coupled chains model introduced by Baeriswyl and Maki ||20|, in 
which, each chain is described by the Su-Schrieffer-Heeger Hamiltonian (SSH) pl| 



H, = -tE[l - (-l)"^,](4+ic,n + H.c.) , (1) 

n 

where j = 1,2 denotes the chain index, and is coupled by an interchain hopping term 

H± = -ts. E(cLc2n + H.c.) . (2) 

n 

Here Cjn is the annihilation operator of the electron at site n on the j'th chain, the spin indices 
have been omitted for simplicity. Each chain is assumed to be dimerized in accordance with 
the Peierls theorem [Q, zj is the dimerization amplitude of the jth chain, \zi\ = \z2\ = z, 
but they may differ in sign. Strictly speaking, the SSH model is directly applicable only to 
polyacetylene, however, recent works have shown that the primary excitation in luminescent 



polymers like PPV can also be described within the linear chain model |123| , |10| , |11| . In PPV 



and its derivatives, the lowest excitonic wave function extends over several repeat units P,p!3 



the properties of exciton are therefore not very sensitive to the delicate structure within the 
unit. From the viewpoint of renormalization [^], we can map the complex structure of PPV 
into an effective SSH system with the same significant physical properties by integrating 
out the freedom of benzene rings and only considering the electrons on the nonbenzene 
carbon atoms [|1^. Thus the features of exciton in SSH model must have some implication 



to that in luminescent polymers. We have also adopted the rigid lattice approximation, 
since many experiments and theories had demonstrated that the primary excitation is the 
exciton and electron-electron interactions are predominant over electron-lattice interactions 



in luminescent polymers [p|,|7|,|8|,p!0|JTl|,|23[] . Theoretical works that have taken into account the 



electron-lattice interaction also show that incorporation of lattice relaxation effect would not 
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lead to an increase in binding energy of exciton . This simplification enables us to handle 
electron-electron interactions in long chains and arrive at an understanding of electronic 
states in luminescent polymers without loss of essential physics, although the quantitative 
explanation of some lattice property like vibronic structure or bond length should, indeed, 
take into account the lattice relaxation effects pB|. 

By introducing operators ajk and hjk through the relation (we take the lattice constant 
a = 1 in this section) |Q 

here is the number of sites per chain, and by using the Bogoliubov transformation 



(3) 



(a \ 




cos 6jk 


sin Ojk ^ 




V J 




^ - sin 9jk 


cos 9jk ^ 


^I3jk J 



(4) 



Hj is diagonalized when tan 26 = Zj tan k 



(5) 



with 



Ek = 2t\Jcos'^ k + Zj sin^ k . 



(6) 



The interchain hopping term then becomes 
H± = -t± Yl[cos{9ik - 02k)ial^a2k + PlkP2k) + sm{9ik - 92k){PlkOi2k - c^lkp2k) + H.c] . 



(7) 



There are two kinds of orderings between the two chains, the parallel ordering {6ik = 62k) 
and anti-parallel one (6*1^ = —62k) HO]- For the case of 9ik = 02k, 



H± = X!("ifc"2fc + PlkP2k + H.c. 



(8) 



the full Hamiltonian H = Hi + H2 + H± can be written as 
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(4k «L) 



2k) 



\^0i2k I 



\ 



-Ek —i± 
-t± —Ek 



1 1 Pik 



y P2k J 



and is readily diagonalized by the orthogonal transformation 



(9) 



{Mk Mk Bik B2kY = Q{0ilk Oi2k Plk (^2kY > 



(10) 



H = YXi^k - - Bl^B^k) + {Ek + t^){A^A2k - B\^B2k)] , (11) 

k 

where, AI^. and Bjj^ {i — 1,2) create an electron in the ith conduction and valence band 
respectively. 

For the case of 9ik = —(^2k = ^k, 

H± = ^[cos 2ek{a\ka2k + PlkP2k) + sin 2^fc(/?|,a2/c - a\kP2k) + H.c] , (12) 



the total Hamiltonian reads 



k 



Ek —t±cos29k i_Lsin2^fc 
— t_L cos 26k Ek —t± sin 26k 

—t±sm26k —Ek —t±cos26k 
t±sm2dk — i_LC0s2^fc —Ek 



\ ( \ 

aik 



tt2fc 
Plk 

y P2k J 



(13) 



Making the orthogonal transformation 0, we obtain the diagonalized Hamiltonian, 

H = Y.[^ik{A\kAik - BlkBik) + e2k{AlkA2k - BlkB2k)] , (14) 



with 



Elk = J El + tl+4:tt± COS k , 



(15) 



£2fe = \jEl + t\- 4tt_L COS k , 



(16) 



and, meanwhile, the transformation matrix 0- 



We add the long range Coulomb interaction He-e to H and study the exciton state, 

He-e = ^ E E yU'PllspJl's' , (17) 
^ IJW ss' 

where pus = c\igCiis — 1/2, I and J are the chain indices, and the interaction potential 

yiJ _ 



Viv , I = J 



The intrachain Coulomb interaction is the commonly used form 

Here U is the on-site Hubbard repulsion, V the nearest- neighbor Coulomb interaction. We 
assume that the interaction potential between the two chains has a relatively simple form. 



iRi+(i- ir- 

Rq will be set to be 1 in following calculations. 

The procedure to determine the exciton state is a standard one. First, we use the single 
particle state of H to construct the ground state l^f), l^f) = Hjk ^jki^jkil^)- Within the SCI 
approximation, the exciton state can be achieved by diagonalizing the Hamiltonian H + Hg-e 
in the subspace of the single electron-hole excitation 

|2, h-j, K) = -L(4,^^5,,„T ± Al^^B^kAM , (18) 

where + is for the spin singlet, — for one of the triplet, kc and —k^ are momenta of the 
electron in conduction band and the hole in valence band, respectively, i and j (=1, 2) are 
the band indices. 

For the spin singlet exciton, 

{i', k'^;j', k'jH + He-e - Eo\i, k^J, K) 
= Sk'^,kJk'^,kASjf{6ii'eik, + E y"{A'kA){cji'c\i){cnAl,^J) 

IJIV 

+ 5u\-5,,>e,k. - E Viy{B],,^cjv){cnc\,){ciB^,^))] + 2Ex - Ec , (19) 

IJW 



where Eq = {g\H + He-e\g), and for the triplet, 

{i', k'^;j', k'jH + H^^^ - £'o|i, k^\j, k,,) 

IJW 

+ 5u\-6,,,e,u. - l^[/(4,„cj,)(c«4,)(cl,5,fc„))] - Ec , (20) 



IJIV 



and 



IJW 



Ec=Y. V!/{A,,.cl,) {cjvAl){B],,,cn) , (22) 

IJW 

where (...) = {g\..\g). We show how to evaluate (cl^AjytJ, (Bj/^^cji), and {c\i,cn) in the 
Appendix. 

This two-coupled-chain system also has the symmetry with respect to the spatial inver- 
sion at a bond center like a single chain. The inversion operator R is defined by 

-R[Cjn] = CjN-n+1 i (23) 

and it is easy to prove 



so 



c^kf^^k. \9) ^ -e^('==-'-")aL,,/,.fe„ 1^7) , (24) 



and from the transformation matrix 0, which diagonalizes the Hamiltonians (P) and (|13|), 
we also have 

^ci j'j k^) > — e " |z, —kc] J, ~ky) . (25) 

Thus we can construct the symmetric state {A) and anti-symmetric one {B). The A 
state is written as 
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\hkciitky] ) — kc', j, ky) " "•'li, kc]j, ky)) , (26) 

and the B state is 

\hkc] 3)ky] +) = —j={\i,kc, i,ky) + e^^^'' ^''^\i,—kc]i,—ky)). (27) 

In the numerical calculation, we may confine ourselves in the A or i? subspaces to diagonalize 
the Hamiltonian H + He-e since the matrix element between A and B states vanish. In the 
exciton state, the relative and center-of-mass motion can also be separated by introducing 
the variable k and K so that kc = k + K and ky = k — K. k and 2K are the momenta of 
the relative motion and center-of-mass motion of the electron-hole pair. Since we are only 
interested in the relative motion, we will just consider the case of K = 0, and now the basis 
is \i,k;j,k;±). 

The wave function of exciton in real space can be determined by 

^l:{I,n■J,l■±)=Y.{I,n,]■J,l,l\^,k■,J,k■±){^,k■J,k■±\^P) , (28) 

ij,k 

the positions of electron and hole are at site n on the Jth chain and at site / on the Jth 
chain, respectively, {i, k; j, k; can be obtained by diagonalizing the matrices (|T9[) and 
(^), and 

(/,n, t; \i,k;j,k;±) 
= ± {cinAU){cl,B,_,)) . (29) 



III. NUMERICAL RESULTS AND DISCUSSIONS 

We have carried out numerical calculations to show the properties of excitons in coupled 
chains. The intrachain parameters are fixed with realistic values, U = 2t, V = t, and 
z = 0.2. The interchain coupling t± can be estimated according to the results of the LDF 
calculations for PPV by Vogl and Campbell In PPV, the total interchain coupling for 
a pair of monomers tmon = 0.64 eV p8|. Naturally, we set the total coupling for a pair of 
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unit in the SSH model 2t± to be tmon- Thus t± for PPV ranges from 0.12t to 0.1 5t. In PPV, 
the nearest interatom distance on adjacent chains is 2.495 A |^8[, while the bond lengths 
are around 1.4 A, so the Coulomb interaction strength U = 0.5t is a reasonable estimate 
for PPV. We will also vary the interchain hopping t± and interchain Coulomb interaction 
parameter U to make the interchain effects more transparent. The system we study consists 
of two chains of = 200. The exciton wave functions in the following figures represent the 
relative distribution of the hole in different positions when the electron is located at site 
n = 101 on the first chain. The 200 + ith site in figures means the ith site of the second 
chain. 

First we study the case of parallel ordering. Figure 1 illustrates the wave functions of 
the lowest and states with t± = 0.15t and If = 0.5t. Since the actual wave function 
contains rapid staggered oscillation with a period of 4a, we have plotted ■?/'(/ = l,n = 
101; J, /; ±)//(/) with /(/) = V2cos[f (/ + ^)] for the odd number / ]18[. We give in Tables 



I and II the exciton energy E^^ and the probability P that the hole is at the second chain 
for several groups of t± and U. /^From these tables, we can see that the hole has more 
probability to stay at the second chain for both B and A states as the interchain hopping t± 
and interchain Coulomb interaction U increase. Since the energies of the A states are close 
to the continuum band and the wave functions nearly extend over the whole system, the 
amplitude of wave function in the second chain is comparable to that of the first chain. The 
energies of B states sit, however, deep in the gap, so that the amplitude of wave function in 
the second chain is smaller and more sensitively depends on the interchain couplings than 
the case of A states. The triplet has similar features of the singlet. 

Now we divert our attention to the case of anti-parallel ordering, which is more likely to 
be realized in practical polymers ||29|. Meanwhile the features of the exciton in this case are 



more interesting. We describe the lowest ^B and -^A states for t± = O.lSt and U = 0.5t in 
Fig. 2. We have plotted ip{I = l,n = 101; J,l; ±)/ f{l) for the odd number / when J = 1 
but for the even number / when J = 2. The exciton energy E^x, the probability P that the 
hole stays in the second chain, and the separation d between positions of maxima of wave 

10 



function amplitude in two chains are listed in Tables III and IV for several groups of the 
interchain parameters. It is shown that the larger U leads to the more chance to find the hole 
in the second chain. However, the dependence of probability P on the interchain coupling 
t± is not so simple as the case of parallel ordering. For the state, P becomes larger 
when increases, while for the state, in contradiction with intuition, P decreases as t± 
increases. Unlike the parallel ordering, where the profiles of wave functions in two chains 
are alike and the maxima of wave function amplitude in two chains sit at the same position, 
in the case of anti-parallel ordering, the wave functions in two chain are well different from 
each other. For states, when the hole approaches the position of electron n = 101 on 
the first chain, the wave function reaches the maximum; but on the second chain, the wave 
function goes to zero when the hole is nearest to the electron. While, for states, when the 
hole and electron overlap on the first chain, the wave function vanishes; but when the hole 
is on the second chain, the wave function reaches the maximum as the hole approaches the 
position of electron. The positions with the largest wave function amplitudes in two chains 
have a separation ~ 5 — 9a for both B and A states. The interchain Coulomb interaction 
can reduce this separation to a certain extent. This big difference between the two ordering 
case is understandable. /^From Eqs. and (|T3|), we can see that, in the case of parallel 
ordering, there is no mixing between valence and conduction bands of different chains and 
the system is more like a two-independent-chain system; in the case of anti-parallel ordering, 
however, the valence- and conduction- band states of these two chains are mixed together 
and the two chains are really coupled. 

For PPV, there are two kinds of ordering between adjacent chains, namely, in-phase and 



out-of-phase ordering [p0|-|32[|. When the practical structure of PPV is mapped into the 
simplified SSH model as stated before, each type of ordering is neither parallel nor anti- 
parallel and the wave function of exciton in PPV should therefore be the composition of 
that in parallel and anti-parallel ordering. For the ^B state, from Figs. 1(a) and 2(a), the 
maximum of the composed wave function in the first chain should sit at the positon of the 
electron, since the profiles of wave function in parallel and anti-parallel ordering are same; 
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while, in the second chain, the position of maximum for the composed wave function must 
be situated between the site facing the electron and the site where the maximum in the 
second chain occurs for anti-parallel ordering. For the state, from Figs. 1(b) and 2(b), 
the composed wave function in the first chain will vanish at the position of the electron, 
and reach the maximum at the site which lies between the position of the electron and 
that of the maximum for parallel ordering. The shape of composed wave function, which 
shows that the position of maximum in the second chain deviates from the location of the 
electron in the first chain, implies that if the interchain exciton is produced, the electron 
and hole tend to be separated with several lattice constant. This is similar as the concept of 
"spatially indirect exciton" proposed by Yan and co-workers to interpret the PA spectrum 
in PPV We emphasize that it is the exchange effect which prevents the electron 

and hole in different chains from approaching each other in the presence of the interchain 
electron-electron interaction, and previous treatments to the interchain Coulomb interaction. 



in which, only the electrostatic energy is included [^], cannot predict this feature. Another 
interesting property is that, for the A state, the interchain exciton is even more likely to 
be created than the intrachain exciton. Recently, it was documented that 80% - 90% 
photoexcitations in PPV are interchain excitations 0, since in practical materials, defects, 
interfaces, and thermal fluctuation can lead to the charge transfer and the mixing between 
the A and B states, the A states with large possibilities of interchain excitons seem to have 
contributions to the great quantity of the interchain excitations in PPV. Obviously, these 
interchain excitations are also important to the PC, which is thought an interchain process 
in PPV. 

In summary, we have studied the excitons in two coupled polymer chains. The wave 
functions of the lowest A and B states for both parallel and anti-parallel ordering have been 
illustrated. We have shown the pronounced difference of the property of exciton between 
these two kinds of ordering. For the realistic PPV, whose properties should be the com- 
bination of that in each ordering, the electron and hole tend to be separated with several 
lattice constant when they are not in the same chain, and in the A state, the probability that 
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the electron and hole are in different chains is even larger than that they are in the same 
chain. These features are helpful to clarify the controversies of interpreting experiments of 
luminescent polymers like PPV. 
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APPENDIX: 

In this appendix, we show how to calculate {cliAi^J, {Bji^^cji), and {cji,cji). It is noted 
that 

{ciiAlj = {n\AljO) , (Al) 



(A2) 



and 



{c]i,cn) = Y.{0\B,,Jjn{Il\BljO) . 



(A3) 



Since from Eqs. (0) and (^), we have 



{2n\aU0) (2n + l|«|j0) 



{2n\(3U0) {2n + l\pU0) 



\ 



( 



Ak2n 



N 



exp{—i9jk) — exp{i9jk + ik) 
i exp{—i6jk) i exp{i6jk + ik) 



(A4) 



with 
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cos^,-, = i=(l + -|)^ , (A5) 

. Zjsink 1 efc i 

here €k — 2tcosk, then combining Eq. (A4) with the transformation matrix 0, we can 
obtain (c|;Afec), {Bjk^cii), and {c^ji,cn). 
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FIGURES 

FIG. 1. Wave functions of the lowest and exciton states for the interchain hopping 
t± = 0.15i and interchain Coulomb interaction U = 0.5t in the case of parallel ordering between 
the two chains, (a) is for state and (b) is for one. 

FIG. 2. Wave functions of the lowest ^B and ^A exciton states for t± = 0.15t and U = OM in 
the case of anti-parallel ordering, (a) is for ^B state and (b) is for ^A one. 
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TABLES 



TABLE I. Interchain parameters {t± and U) and resulting exciton properties (£^ex and P) of 
the lowest states for the case of parallel ordering between the two chains. 



I±. U (I) 


0.03, 0.5 


0.03. 1.0 


0.09, 0.5 


0.09. 1.0 


0.15, 0.5 


0.15. 1.0 




1.095 


1.057 


1.010 


0.941 


0.902 


0.821 


P 


0.062 


0.300 


0.241 


0.428 


0.330 


0.457 



TABLE II. Interchain parameters {t^ and U) and resulting exciton properties (-Eex and P) of 
the lowest states for the case of parallel ordering between the two chains. 

tx, U (t) 0.03, 0.5 0.03, 1.0 0.09, 0.5 0.09, 1.0 0.15, 0.5 0.15, 1.0 
E^^ (t) 1.345 1.331 1.229 1.211 1.109 1.091 

P 0.335 0.163 0.146 0.488 0.168 0.493 



TABLE III. Interchain parameters {t± and U) and resulting exciton properties (-Eex, P, and 
d) of the lowest states for the case of anti-parallel ordering. 

t^, U {t) 0.03, 0.5 0.03, 1.0 0.09, 0.5 0.09, 1.0 0.15, 0.5 0.15, 1.0 
E^ex (t) 1.110 1.096 1.110 1.095 1.108 1.093 

P 0.005 0.007 0.045 0.059 0.111 0.139 

d {a) 7 7 7 7 7 7 



TABLE IV. Interchain parameters {t± and U) and resulting exciton properties (£^ex, P, and 
d) of the lowest states for the case of anti-parallel ordering. 



i.±. u (/) 


0.03. 0.5 


0.03, 1.0 


0.09. 0.5 


0.09, 1.0 


0.15. 0.5 


0.15, 1.0 


£^ex {t) 


1.345 


1.194 


1.328 


1.188 


1.308 


1.178 


P 


0.913 


0.987 


0.701 


0.903 


0.609 


0.798 


d (a) 


9 


7 


9 


7 


7 


5 
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